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We map the problem of self-avoiding random walks in a O solvent with a chemical potential for 
writhe to the three-dimensional symmetric f/(A r )-Chern-Simons theory as N — * 0. We find a new 
scaling regime of topologically constrained polymers, with critical exponents that depend on the 
chemical potential for writhe, which gives way to a fluctuation-induced first-order transition. 
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The statistics and properties of random walks are central to our understanding of a large class of phenomena in 
mathematics, biology, physics and even economics. The broad applicability of random walks in modelling physical 
situations is, in part, due to the Gaussian statistics of its correlations, rendering it attractive from an analytical 
viewpoint. In polymer physics, self-avoiding random walks model polymers in a good solvent that are much longer than 
their persistence length. This identification establishes a theoretical basis for the computation of scaling exponents 
which are, in turn, related to the exponent v that relates the radius of gyration Rq to the length of the chain L 
through Rq ~ L v . In this letter we study the scaling behavior of closed polymers with a chemical potential for their 
writhe Wr. Closed loops of DNA, as found in bacterial plasmids, experience a topological linking number constraint 
enforced by the twist rigidity of the double backbone. Refinements can be made to de Gennes' mapping to incorporate 
this additional physics through the addition of a minimal coupling to a Chern-Simons gauge field [1]. Specifically, 
the N — ► limit of a U(N) Chern-Simons theory with the gauge field coupled to the common U(l) can be used to 
study self- avoiding walks with this topological constraint. The Abelian Chern-Simons gauge field is used to enforce 
the topological constraint in the form of White and Fuller's [2, 3] famous relation [4, 5, 6, 7] Lk = Tw + Wr. It 
relates the topological linking number Lk to the total amount of twist Tw and writhe Wr found in, for example, DNA 
backbones. We find that the topological constraint leads to a new universality class of random walks. 

Without a topological constraint, the individual monomers suffer two interactions: a chemical interaction which 
depends on the solvent and a hard-core interaction which is entropic. In the case of good solvents, these interactions 
of chain connectivity and excluded volume are described by de Gennes' mapping [8] of the 0(A)-symmetric 4 -theory 
in the N — ► limit. The Wilson-Fisher fixed point controls the scaling in d = 4 — e dimensions where v as 0.588 in 
d = 3 [1]. In poor solvents, the polymers collapse into compact globules and Rq ~ . The transition between the 
two solvent regimes is the point which is described by an 0(A r )-symmetric 6 -model at the point where the </> 4 
term vanishes or, in other words, at the tricritical point [9]. Logarithmic corrections to mean-field results have been 
calculated, but it has been proven difficult to test these predictions experimentally [10, 11]. Computer simulations 
[12, 13, 14] have also been performed but agreement with theory has only been tentative. 

The topological constraint imposed by Lk is inherently three-dimensional as is the Chern-Simons term which 
surrogates for it. As a result, the perturbative study of a polymer in a good solvent cannot be performed via a 
conventional e = 4 — d expansion. In [1] it was shown that the exponents at the d = 3 Wilson-Fisher fixed point were 
unchanged at one loop order. Because there is no controlled approximation to go to higher order, we have focussed 
here on polymers at the point with a topological constraint so that a systematic expansion can be performed directly 
in d — 3 — e dimensions. 

At the point, the exponent v is the correlation length exponent at the tricritical point [i = u = of the following 
free energy density for the A-component, complex scalar <f> [8, 9]: 



in the N — > limit. In three dimensions, renormalization group analysis shows that although there are logarithmic 
corrections to mean-field behavior, the scaling exponents agree simply with predictions from dimensional analysis 
[15, 16]. To incorporate the topological constraints found in closed polymers, we introduce an Abelian Chern-Simons 
gauge field with the free energy density [1]: 
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Since we are interested in the critical behavior of our theory near the point, we perform our analysis at the tricritical 
point where /i and u vanish. The energy is the sum F — Fes + Fy } + F g f where 



Fes = \ I ^xe^pA^Ap (3) 
F h = [ d^ld^i-igoAp^f + V^i) (4) 



F 



gf = d 3 x(dM 2 (5) 

V(tfH) = Mo|0| 2 + wo (H 2 ) 3 (6) 

We impose the Landau gauge (A — > 0) in F g f for all our subsequent calculations. Because the graphs first diverge at 
two- loops, standard dimensional regularization is adequate and leads to - poles in d = 3 — e dimensions. To evaluate 
each two-loop diagram, we perform the first momentum integral in d = 3 and then use dimensional regularization 
on the remaining single integral. This scheme works because the first integral can only give power law divergences 
since any integrand with odd powers of momentum vanishes, leaving only even powers of the momenta compared to 
the three-dimensional measure. Moreover, by power-counting, the remaining integral must diverge logarithmically 
and cannot cancel the power-law divergence from the first integration. Thus the logarithms only arise in the second 
of the two integrations. A more sophisticated treatment suggests that this scheme is consistent to all orders [17]. 
As is usual, we first perform the necessary tensor algebra in physical dimensions before analytically continuing the 
dimensions of the resulting scalar integrand. The correlation functions satisfy the Slavnov- Taylor identities [18] and 
the Ward identities [17]; gauge invariance is therefore preserved. 

We first introduce the renormalized free energy as follows: 

+ l -Z A e^Apd v A p + ^{d^f + Z w w (V| 2 ) 3 }, (7) 

where g = gM^ 3 ~ d ^ 2 , fl = /iM 2 , w = wM 6 ~ 2d , and M is the momentum scale at which we are renormalizing. 
Denoting Y^ 
equation is: 



Denoting Y^ 4 "^^ as the renormalized proper vertex of scalar fields and Na gauge fields, the Callan-Symanzik 



where M is the renormalization scale. The engineering dimension of any correlation function Y depends on its field 
content so that d(Y) = 3 - ±iV - Na- We have 

m = M ^ lnZ <P ( 9 ) 

VA = M^\nZ A (10) 

Mo = M M 2 §i (11) 

50 = gM ^-^-^j= (12) 

w = wM 6 ~ 2d ^ (13) 



Perturbatively, the underlying symmetry of our action provides a tight constraint on the divergences in our calcu- 
lations. The odd parity of the Chern-Simons field prevents any correction to the scalar field at the one-loop level: 
by rescaling A^ — > g~ 1 A fll we see that under parity g 2 — ► — g 2 , and so correlations of <f> which are parity invariant 
can only depend on g 4 and thus the first corrections are at two loops. Accordingly, we do not find any one-loop 
contributions to or Z^. Further, the Coleman-Hill theorem shows that the j3- function of the Chern-Simons gauge 
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FIG. 1: Two-loop diagrams arising from the gauge field coupling to the scalar field: contributions to (a),(b); corrections to 
the cubic gauge vertex (c)-(f); contributions to Z M (p=0) (g). 



TABLE I: Divergent contribution from each graph and to the appropriate renormalization constant. Note that there are no 
divergences at one loop. 
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coupling receives no contribution beyond one-loop [19, 20, 21, 22] in perturbation. More simplification results from 
the vanishing of diagrams at any order with closed scalar loops: they necessarily introduce a combinatoric factor of 
N and do not provide corrections to the gauge field as TV — > and so a Maxwell term, F 2 though allowed by sym- 
metry, is not generated. The details of the combinatoric argument go as follows: consider a graph with only external 
gauge-field legs. Since an external gauge-field leg must necessarily connect to two internal 4> legs (and possibly one 
internal gauge-field leg) and since the U (N) index of is not carried by any of the other external legs, there must be 
a sum over that index. Since that sum is proportional to N, this graph vanishes as N — > 0. More complex internal 
topologies only add factors of N to the graph and do not change this result, thus Za = 1. For the same reason there 
exists only a handful of potential two-loop contributions to the renormalization functions Zx ■ 

All the non-vanishing two-loop Feynman diagrams are shown in Figures 1 and 2. Figs, la and lb are contributions 
to Z^, while Fig. lc through If are corrections to the cubic gauge vertex. Fig. lg is the same as Fig. la but 
evaluated at zero external momentum, contributing to Z^. The evaluation of these graphs is straightforward, since 
their singularity structures involve only simple poles by power counting. Our results for each graph agree with those 
in [17]. 

Employing the results in the table, we have 
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FIG. 2: Two- loop contributions to r' 6 '°\ Figures (c), (d) and (e) are representative of the topology; we have accounted for 
other contractions of 6 with d>* . 



from which we find that at two-loop order in d 
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where 
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For nonvanishing 5, as e — > both roots are complex and there is no physical w fixed point. However, for g 4 < 0.851e, 
the radicand in (23) is real and w + is a stable fixed point for e > 0. Note, however, that for e > the gauge coupling 
runs away to large values and thus the only stable fixed point in d = 3 — e is at (g, w) — (0, 47r 2 e/33). Focussing on 
d = 3, we see that g is exactly marginal and w± = — g 4 (0.242 =p 0.392i). Writing M — M e~ l , we have 



w(£) = 3?w + + 3w + cot ( cot 



3w + 



(24) 



where n 



2^ 



33 



= 0.655. Note that this solution is unstable and as I grows runs away to negative values of 



10, reminiscent of the behavior in a superconductor [23] and signaling a first-order transition. However, until any 
new fixed point controls the scaling, the critical exponents v and 77 are determined entirely by the exactly marginal 
coupling g: 
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We thus see that a chemical potential for writhe can alter the radius of gyration exponent v and therefore writhe 
alters the universality class of a self-avoiding walk before driving it to collapse. The gauge field is not perturbatively 
renormalized and thereby preserves its topological character. As in [1] the scaling behavior of the average writhe (Wr) 

and the average squared writhe (Wr 2 ) are given in terms of the specific heat exponent a = 2 
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we see that if the chemical potential vanishes then (Wr) — as expected. These logarithmic corrections are in 
addition to the writhe that is stored is polymer segments comparable to the persistence length (e.g. plectonemes) 
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which scales as L, a result which does not violate the rigorous bound ln(|Wr|) > 5I11L [24, 25]. Our expressions 
of scaling exponents rj and v depend continuously on g 2 and are thus similar to those in the two-dimensional XY 
model. There, topological defects are responsible for this uncommon behavior; here the topological link constraint is 
responsible. 

In conclusion, we have mapped the study of the statistics of self-avoiding random walks in a 8 solvent to a Chern- 
Simons field theory and have found a new scaling regime for topologically constrained polymers by calculating the 
scaling exponents r\ and v to two-loop order at the 6 point. We conjecture that even in a good solvent the scaling 
behavior will be altered, though as prior work indicates [1] it is difficult to establish results in a controlled approxi- 
mation. Since taking N — > amounts to canceling a functional determinant, progress in this problem might be made 
by introducing fermionic partners to the complex scalars to have the same effect. It is possible that supersymmetric 
formulation of this field theory would yield a more complete understanding of the effect of topological constraints. 
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